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We propose that flow fluctuations have the same origin as transverse momentum fluctuations. 
The common source of these fluctuations is the spatially inhomogeneous initial state that drives 
hydrodynamic flow. Longitudinal correlations from an early Glasma stage followed by hydrodynamic 
flow quantitatively account for many features of multiplicity and pt fluctuation data. We develop 
a framework for studying flow and its fluctuations in this picture. We then compute elliptic and 
triangular flow fluctuations, and study their connections to the ridge. 
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I. INTRODUCTION 

Fluctuations in the early stages of nuclear collisions 
contribute to the anisotropic flow measured in RHIC and 
LHC experiments, particularly the odd harmonics. In 
Ref. pQ we find that this early-time variation can also 
account for the multiplicity and transverse momentum 
fluctuations measured in the same experiments. Never- 
theless, flow and fluctuation observables reveal different 
aspects of the initial state. Harmonic flow is generated 
largely by the global spatial anisotropy of the initial state. 
In contrast, fluctuation observables probe local spatial 
correlations that are largely independent of the overall 
geometry pQ. 

In this paper we argue that the fluctuations of the har- 
monic flow coefficients are driven by the very same local 
correlations that give rise to pt and multiplicity fluctu- 
ations. Moreover, if early time correlations are the only 
source of these fluctuations, then the two measurements 
are related with no free parameters. 

Experimenters define flow fluctuations by studying the 
difference of the flow coefficients w„{2} and w„{4} mea- 
sured using two and four particle correlations, respec- 
tively Jets, resonance decays, and HBT effects also 
contribute to that difference. Measurements using parti- 
cles separated in rapidity by more than 1 — 2 units elim- 
inate these short range effects. Significantly, a difference 
between v n {2} and v„{4} remains [BJ, suggesting that 
flow fluctuations are a long range phenomenon. Causal- 
ity dictates that such long range correlations originate in 
the early stages of the collision [7j |8] . 

To compute flow fluctuations we build on an approach 
started in Refs. [HI H] in which long range correlations 
result from the fragmentation of Glasma flux tubes. In 
this formulation local spatial correlations emerge from 
fluctuations in the number and distribution of flux tubes. 
These spatial correlations are then modified by transverse 
expansion, giving rise to azimuthal correlations. A simi- 
lar physical picture motivates studies using a wide range 
of different techniques [T0MT7] . 

This paper is organized as follows. In SecjllJ we be- 
gin with a general description of multi-particle correla- 



tions and flow coefficients. We relate the intrinsic correla- 
tions in the multi-particle system to flow and its fluctua- 
tions. To begin, we write expressions for v n {2} and v n {&} 
in terms of correlation functions, adapting the tools in- 
vented in Refs. |T51 [TS] to our framework. Further results 
are in the appendix. We next define flow fluctuations in 
terms of v n {2} and u n {4}, and turn to discuss the phys- 
ical interpretation of this definition. 

We discuss the common influence of local correlations 
on flow and pt fluctuations in Sec |III| This relationship is 
the heart of our work. The fluctuations of these quanti- 
ties are both consequences the spatially inhomogeneous 
collision environment modified by hydrodynamic flow. 
The general arguments in Sees. [Il] and |III| set the stage 
for the more phenomenological analysis that follows. 

In Sec lIVI we describe initial state fluctuations in a 
CGC-Glasma picture. The number and distribution of 
Glasma flux tubes relative to the geometrical shape of the 
system ultimately determine the energy, projectile-mass, 
and centrality dependence of flow fluctuations. Next, in 
SecjVJwe explain how collective flow and the correlation 
to the reaction plane modify local correlations. 

We calculate the contribution of long range correla- 
tions to the flow coefficients and their fluctuations in 
SecfV] Our results are in good agreement with the latest 
LHC and RHIC data. In Sec |VI| we argue that these 
same flow fluctuations are also responsible for the ridge. 
The same factors that influence flow fluctuations also 
determine the multiplicity and momentum fluctuations 
computed in pQ. We emphasize that no new parame- 
ters are introduced here, so that results can be compared 
directly with pQ. 

In Sec I VIII we discuss the extent to which the flow co- 
efficients determine the angular distribution of the ridge. 
We also discuss the possible factorization of the Fourier 
coefficients of the pair distribution into products of flow 
coefficients [2"Ull23| . We show that flow fluctuations can 
violate factorization at low p t , depending on how the flow 
coefficients are defined. We then speculate that factor- 
ization holds for momenta > 1 — 2 GeV regardless of the 
source of anisotropy. 

Lastly, we summarize and discuss the broader implica- 



2 



tions of our results in Sec lVIIll 



II. FLOW AND ITS FLUCTUATIONS 

Nuclear collisions at non-zero impact parameter b pro- 
duce anisotropic flow [SH 123- This anisotropy derives 
from the change in the shape of the collision volume with 
respect to the reaction plane, i.e., the plane spanned by 
b and the beam direction. If the reaction plane is known, 
this anisotropy is characterized by the moments 



(cosn(0- ip RP )}, 



(1) 



where is the azimuthal angle, ip RP is the angle of the 
reaction plane, and the brackets denote an average over 
particles and events. While the reaction plane cannot 
be observed directly, its influence can be deduced from 
multi-particle correlation measurements. Many strate- 
gies have been employed for measuring these flow coeffi- 
cients, and they all have strengths and weaknesses. 

Experimenters often deduce the flow from the two- 
particle cumulant 



v n {2} 2 = (cosn(0i - 2 )), 



(2) 



exploiting the reaction plane information implicit in the 
relative distribution of particle pairs. Specifically, they 
measure the relative azimuthal angle for each particle 
pair in each event, and then average over events to obtain 



,{2} 2 = 



(E«,cosn(^, - <f>j)) 



(N(N 



1)> 



(3) 



|26j . The cumulant method was developed in Ref. [T51[T§] 
and has seen extensive use [211 ES] . 

In this section we cast the framework of [T51 Q15] in 
terms of our formulation of correlations in pQ to study 
event-wise flow fluctuations. Many results in this section 
are implicit in |18| 119] ; only our application to the com- 
putation of fluctuations is new. We start by writing ([3]) 



V n {2} 2 



P2(Pl,P2) 

(N(N-l)) 



b 2 )dp 1 dp 2 , (4) 



where the integrals are over the momenta of particles 
for i = 1, 2 and the distribution of particle pairs is 



P2(Pl,P2) 



dN 



dpidp 2 



(5) 



In the absence of correlations, P2(Pi,P2) — > 
/°i(Pi)Pi(P2)i where the single particle distribution 
is pi(p) = dN/dp. We stress that the densities p\ 
and p 2 are event-averaged quantities that respect the 
reaction plane. The factorization of p 2 then allows for 
factorization of Q such that v n {2} 2 — > (v n ) 2 . 

Including two-particle correlations, the pair distribu- 
tion does not factorize. To identify the contributions of 



the mean anisotropic flow ([I) and genuine two-particle 
correlations to v n {2}, we fo low [18] and write ([5| in 
terms of a cumulant expansion, 



P2(pi,P2) = Pi(pi)pi(P2)+r(pi,p 2 ), 



(6) 



where r is the two-particle correlation function. We then 
use 0, Q, and ([6} to write 



v n {2} 2 = 



2at 



(7) 



This result is standard, and the factor of two in a is 
conventional [27] . The contribution to ^ from two-body 
correlations is 



v n = / dp!dp 2 



r(Pi,P2) 

2(N(N -1)) 



cos nA0, 



(8) 



where A0 = 0i — 02 is the relative azimuthal angle of 
the correlated particles. 

Strictly speaking, this quantity measures the event- 
wise fluctuations of anisotropic flow due to geometry and 
dynamics together with "non-flow" correlations from res- 
onance decays, the HBT effect, and jets. We will not 
address non-flow effects here. Observe that all such cor- 
relation effects are of order 1/N, where N is the number 
of particles in the system. We focus on contributions to 
fluctuations at leading order in N. Accordingly, we omit 
a factor (N) 2 /(N(N - 1)) « 1 that multiplies (v n ) 2 term 
in Q; see appendix |A| for exact formulae. 

To reduce the effect of fluctuations on the flow signal, 
experimenters also measure the four particle cumulant 
w„{4} defined by the relation 



v n {4} 4 = 2v n {2} A - (c 



+ 02-03-04)) (9) 



[TBI [12] . The four-angle term depends of the four particle 
correlation function p 4 (l,2,3,4) = P4(pi, Pa, P3, P4)- In 
systems with large numbers of particles, multi-particle 
correlations are dominated by two-particle correlations. 
We therefore write 



p 4 (l,2,3,4)=p 1 (l)p 1 (2)A> 1 (3)pi(4) 
+ p 1 (l)p 1 (2)r(3,4) + 
+ r(l,2)r(3,4) + ... 



(10) 



where pi(l) = Pi(pi), etc., and the ellipses represents 
the distinct permutations of the momenta. Computing 
(cosn(0! + 02 — 03 — 04)) j one finds that the second and 



third lines in ( 10 1 contribute respectively terms of order 
{vn) 2 &n and (T^~that precisely cancel the a n contributions 
from v n {2} 4 . This cancellation is by design - it dictates 
the form of Q proposed in Ref. [TH1 QI5] . One finds 

Vn{±} W («„), (11) 

plus corrections that are in practice very small. The 
fact that i>„{4} receives no contribution from a n is usu- 
ally thought of as a consequence of the Bessel-Gaussian 
approximation ]27] . but we see that it is true whenever 
two-body correlations are dominant. However, note that 
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FIG. 1. Measured v 2 {2} and v 2 {4} from STAR g] compared 
to calculations. Computed «2{2} (solid) uses and (S), 
while «2{4} (dashed) uses Q, mf, and (24}. The difference 
between the calculated curves is due to flow fluctuations. 



definitions of (v n ) other than ([I]) would yield different 
results. 



As an aside, we briefly discuss the corrections to (11 1 
Strictly speaking, calculation of ([9| using ( 10 1 yields 



«„{4} 2 « {v n f 



2£f 



(12) 



where 



K = J d Pl dp 2 T^j^ cos 2n($ - ) (13) 

for $ = (0i + <j>2)/2 the average angle of the correlated 
pairs. Further corrections of higher order in 1/N are 
discussed in appendix |A") Both S n and <7„ are of the same 
order in 1/N, but one expects E ra to be smaller than a n 
because it effectively is a higher harmonic, of order 2n 
rather than n [X5] . Indeed, we find the contribution of 
£„ to u n {4} smaller than 1.2% for n = 2 for peripheral 
collisions in our model. Note that £„ = in central 
collisions. We comment that E n represents fluctuations 
of harmonics of which are important in studying the 
Chiral Magnetic Effect [25] . 

Our aim is to understand the contribution of two- 
particle correlations to flow fluctuation measurements. If 
the only correlations are due to the reaction plane then 
v n {2} — v n {4}. This requires that no other sources of 
correlation exist, i.e., r(p 1 ,p 2 ) = 0, allowing P2(Pi7P2) 
to factorize. However, experiments measure differences 
between v n {2} and w n {4} for several harmonic orders, 
n [2H5]. For example, Fig [I] shows i>2{2} and «2{4} as 
measured by the STAR experiment in 200 GeV Au+Au 
collisions [5] . One can quantify this difference in terms of 
er„ by substituting the square of (11) into ([7]), to obtain 



v n {2} 2 - v n {AY 



(14) 



Measurements with rapidity separations should suppress 
contributions to er„ from jets, resonance decays, and 



other short range non-flow correlations, but do not ex- 
plain the difference between v n {2} and u„{4}, suggesting 
that (w^)-factorization is not a signature of flow or flow 
fluctuations. 

An alternative measure of flow fluctuations can be con- 
structed using event-by-event harmonics [37J [2H1 130] • If 
each event produces a set of harmonics {v n }, then one 
can compute the variance a 2 , = (vf/) — (v n } 2 for an ensem- 
ble of such events. It is then reasonable to ask how our 
a n compares to this variance. Taking the mean square of 



V n to be (vl) = ((£. 



cosnl 



«V)) 2 )/W 2 , we find 



S 2 /2, 



(15) 



where a n and £„ are given by (8]) and ( 13 1, and we keep 
only leading order in N . We define the statistical vari- 
ance 



1 



[V27-, 



2 (TV) 



(16) 



These fluctuations arise because a finite number of parti- 
cles are sampled in each event. In contrast, the dynamical 
fluctuations a n and £„ are caused by the correlations be- 
tween particles 31]. To an excellent approximation the 
V2n contribution to ( 16 ) and the contribution to ( 15 ) 
can be neglected. 

We see that the event-wise variance of v n is very dif- 



ferent from the quantity ( 14 1 that characterizes the dif- 



ference between v n {2} and w„{4}. The statistical contri- 
bution <7 s tat can be comparable in magnitude to a n . We 
point out that the quantity a n may be measured directly 
from event-by-event harmonic coefficients. 

It is instructive to compare our results to the model of 
Ref. [27l|30j in which flow fluctuations are exclusively due 
to the event-wise variation of the global geometry. In the 
purely geometric interpretation, one defines an eccentric- 
ity for each event e n . If one assumes that the relation be- 
tween e n and the resulting anisotropy of the fluid flow v n 
is approximately deterministic, then fluctuations of the 
ratio v n /e n are negligible. The geometric contribution 



to the variance is then a^/{v n ) 



l )l{Zn) 



We mention that Glauber model estimates of the relevant 
eccentricities are consistent with data for n — 2, where 
the correlation of v n with an event-wise eccentricity is 
highly plausible and well established [32]. However, the 
applicability of such an approach for n ^ 2 is less clear. 

We comment that the purely geometrical description 
of flow fluctuations in Ref. [30] requires (v n ) ^> a n in ad- 
dition to v n oc e n . Neither our general discussion here nor 
our specific calculations in Sec. |V| require these assump- 
tions. In Sec. [V] we include geometric fluctuations in our 
estimates of flow fluctuations through the distribution of 
flux tube sources p FT . Correlations are then modified by 
the local transverse expansion of the system. 
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III. LOCAL CORRELATIONS 



IV. SOURCE OF CORRELATIONS 



To understand the source of flow fluctuations, we recall 
some lessons from general fluctuation studies [3TJ |3"3"H3"5] . 
Such studies typically focus on the variation of bulk ob- 
servables such as multiplicity or transverse momentum 
within an ensemble of collisions. We draw most heavily 
from the study of p t fluctuations in Ref. [T] , to which the 
present work is essentially a sequel. 

We measure pt fluctuations using the covariancc 
(8pti8pt2), where 5p u = pu - (pt) is the deviation of 
each p t from the average. Pairs in which both particles 
have higher than average momentum add to (SptiSpts)- 
Lower-than-average pairs also add to the covariance, 
while high/low pairs subtract from it. In global equi- 
librium (6ptiSpt2) = 0. The presence of hot spots makes 
{8pti8pt2) > (as would cold spots) [35]. Motion of the 
sources further enhances this quantity [TIT]- It follows 
that both jets and flow add to the p t covariance. 

Moving hot spots in the fluctuating system are sources 
of local correlations. We distinguish such correlations 
from those due to the global event-wise variation of the 
shape and size of the collision volume. We write the p t 
covariance as 



(SptiSp t2 ) = / dpidp 2 



KPl,P2) 

(N(N-l)) 



SptlSp: 



12- 



(17) 



In [T] we demonstrate that the p t covariance is indepen- 
dent of the global spatial anisotropy, because (17) is in- 
tegrated over azimuthal angle. It is also independent of 
size fluctuations by construction [37 . 

To see that flow fluctuations are also primarily driven 
by local correlations, observe that ( 17) has the same form 
as ([8| with 5pt\8pt2 replaced by cosnA0. As the pt co- 
variance probes correlations in momentum, so a n mea- 
sures correlations of the angular separation A</> of pairs. 

This analogy makes physical sense because of the com- 
mon influence of hot spots on A(j> and p t . Flow pushes 
particles from a source into a particular opening angle, 
depending on the source position and the local fluid ve- 
locity. The same flow velocity boosts the pt of these 
particles. Therefore multiple random hot spots give rise 
to both flow and p t fluctuations. Based on this analogy, 
we expect a n in central collisions to be non-zero for all 
orders n, while (v n ) defined by ([!]) vanishes. Further- 
more, it is not likely that a n is as strongly influenced by 
the global geometry as the mean (v n ) . 

In the next sections we discuss flow fluctuations due 
to early time local spatial correlations contributing to 
r(pi, P2). We employ a Glasma-flux tube model in which 
particles are initially correlated at the point of produc- 
tion and modified by later stage transverse expansion. 
These flux tubes produce our 'hot spots'. Transverse 
expansion is modeled with a blast wave scenario which 
inherits anisotropy from the average eccentricity (e RP ). 



Flow transforms early-time spatial correlations into 
observable momentum correlations. Two-particle corre- 
lations probe the initial state average pair distribution 
ri2(xi,X2), which encompasses both genuine pair corre- 
lation as well as correlations due to fluctuating event ec- 
centricities, shapes, and orientations (with respect to the 
reaction plane). Define the initial state spatial correla- 
tion function 



c(xi,x 2 ) = n 2 (xi,x 2 ) - ni(xi)m(x 2 ) 



(18) 



where rii(xj) is the particle density. In the absence 
of genuine two-particle correlations the pair distribu- 
tion factors into a product of one body distributions, so 
that the correlation function vanishes. Correlations come 
about because the position and number of the colliding 
nucleons fluctuates in individual events. In general, these 
fluctuations result in high-density 'hot spots' of vary- 
ing size distributed throughout the collision volume. It 
is more likely to find pairs of particles near hot spots. 
Averaging over an ensemble of events therefore yields a 
nonzero c(x!,x 2 ). We define the correlation strength 1Z 
as / c(xi,x 2 )d 3 xid 3 x 2 = (N) 2 TZ. 



In Ref.[I] we construct a correlation function (18) to 
describe a specific early-time scenario based on a Glasma 
theory. In that picture parton production arises from the 
fragmentation of longitudinal color fields created at the 
moment of the collision. It is helpful to think of the 
longitudinal fields as flux tubes of transverse size ~ Qj 1 , 
where Q s is the saturation scale. The saturation scale 
depends on many collision variables including the density 
of participant nucleons and the collision energy [38l [39] . 

Particle production then proceeds via a system of flux 
tubes distributed over the collision volume. Given that 
the transverse area of a flux tube is much smaller than the 
transverse collision area and taking particles emerging 
from the same tube as correlated, we write 



c( Xl ,x 2 ) = (N) 2 KS(r t )p FT (H t ). 



(19) 



The spatial coordinates are relative, r t = r fl — r t2 , and 
average, R t = (r t -y + r t2 )/2, transverse elliptical coordi- 
nates. The delta function enforces the correlation due 
to the common point of production, and p FT (Rt) is the 
average transverse distribution of flux tubes 



p FT (R t ) 



(20) 



Here the shape of (20) resembles the nuclear thickness 
function and ttR\ is the area of the overlap region. The 



distribution ( 20 ) represents an average over all possible 



shapes, which we have taken to be a simple ellipse. 

To compute the correlation strength 1Z, we imagine 
each event produces K flux tubes with transverse size 
~ Qj 2 - In the saturation regime K is proportional to 
the transverse area R\ divided by the area per flux tube, 
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Q s 2 [HE] • Allowing K to fluctuate from event to event 
with average (K), we calculate in Ref.[8] that 



n 



(N 2 ) - (N) 2 - (N) 

W 2 



cx (KY 



(21) 



Each Glasma flux tube yields an average multiplicity of ~ 
&J 1 {Q a ) gluons and as in Ref.[58]. the number of gluons 
in a rapidity interval Ay is then 



(N) = (dN/dy)Ay ~ a.-\Q.)(K). 
Finally, the Glasma correlation scale, 
UdN/dy = Ka- 1 (Q 2 s ), 



(22) 



(23) 



follows from the multiplication of (21 1 and (22 1. Notice 



(23 



sea 



) is dimensionless and depends only on the saturation 
e, Q 2 , which can be calculated from first principles. 
Measurements of the ridge at various beam energies, tar- 
get masses, and centralities fix the dimensionless coeffi- 
cient k and are in excellent accord with the leading-order 
dependence 0[H]. 

In the next section, we will see how transverse expan- 
sion modifies Glasma correlations (19), but it is impor- 



tant to note that ( 23 ) significantly impacts the centrality 



dependence and determines virtually all of the energy 
dependence. 



V. FLOW FLUCTUATIONS FROM GLASMA 
AND TRANSVERSE EXPANSION 

Transverse collective expansion gives rise to momen- 
tum space correlations in two ways. First, the event 
geometry influences global correlations; collective flow 
emerges with respect to a reaction plane. Second, the 
local fluctuations, that together determine the global 
shape, individually induce local correlations; flow, on av- 
erage, modifies the momenta from all partons emerging 
from the same fluctuation in a common way. 

To understand how flow modifies local correlations 
we first model the average anisotropic expansion of the 
system with a blast wave model [H [9l I40H42] . As in 
Refs.0[9] the invariant single-particle distribution, 



pi{p) ^dWirI f ^ p)dT ' 



(24) 



is defined as the particle flux through a surface, a 
[43] . At constant proper time, Tp, the flux element 
is dT = p^da^ = Turrit cosh(y — r])drid 2 r, where 
7] = (l/2)ln((t + z)/{t — z)) is the spatial rapidity. 
Here /(x, p) = (2ir)~ 3 exp{— u^p^/T} is the Boltzmann 
phase-space density for a temperature T and fluid four- 
velocity U/j,. We follow Ref-HQ] and write the four ve- 
locity of the longitudinal-boost invariant blast wave as 
u n = 7t( c °sh?7, v t , sinh?7), with v t , the transverse veloc- 
ity and 7( = (1 — v 2 )~ x l 2 . The phase space density is 
then / cx exp{-7 t m t cosh(j/ - rj)/T} exp{7 t v t • p*/T}. 



To account for the initial elliptical shape and 
anisotropic velocity we require one additional parame- 
ter, the (average) event eccentricity, e. The anisotropic 
transverse flow velocity, 



7 t v t = A (e K x + e y y) 



(25) 



mimics the geometrical parton density gradients. To ac- 
count for the larger pressu re in the direction of the reac- 
tion plane (x- dire ction in (25)), we take s Xj y = yl ± e. 



Rewriting Eq.(25 



as (j t v t 



■r 2 {l + ecos(20)), 



can see how this choice of flow velocity results in an asym- 
metric flow velocity peaked in the direction of the reac- 
tion plane. Note, r is now the elliptical polar radius and 
depends on e and the spatial angle. 

In this formulation, particles passing through the 
freeze out surface have a velocity, v s , that is connected to 
the position of production, thus the transverse velocity 
is constrained via jfv 2 — ^ 2 R 2 ur f ace - As in p], we take 
the blast wave parameters v s and T from experiment |44j 
and parameterize e to fit i>2. However, in Ref. [T], we did 
not include the impact of fluctuations on the calculation 
of w n -like correlations - a focus in this work. Hence, in 
Ref.pQ, e is chosen to fit V2{2}. In this work we identify 
(v2) BW with purely reaction plane angular correlations, 
more suitably comparable to «2{4}. In Figfl] the dashed 
line shows this comparison where the open squares are 
the STAR Au-Au 200 GeV w 2 {4} measurement 0. 

Now that our blast wave formalism is tuned to re- 
produce the average transverse expansion in heavy ion 
collisions, we are in a position to calculate the momen- 
tum space two-particle correlation function. As discussed 
by other authors, see e.g. [45H47] . the two-particle joint 
probability distribution ([5| which we rewrite as 



r(Pl,P 2 ) = P2(Pl,P2) - Pl(Pl)pl(P2) 



(26) 



to highlight the fact that any genuine correlations, in- 
cluding those from local initial state fluctuations (18), 
break the factorization P2(Pi,P2) — > Pi (pi)pi(P2)- 



The spatial correlation function, (19), embodies the 



correspondence between the two particle joint probabil- 
ity distribution and transverse expansion. Intrinsically 
correlated partons originate from the same transverse po- 
sition and experience, on average, the same transverse 



momentum modulation from flow. Generalizing (24), we 
represent the genuine two-particle correlation as 



r(Pi,P2) =j c(x 1 ,x 2 )/(x 1 ,p 1 )/(x2,p 2 )dr 1 dr 2 



(27) 



Notice that full integration of (27) yields (N) 2 1Z as does 
(18), tying azimuthal correlations to multiplicity and 
transverse momentum fluctuations pQ. 

To compute flow fluctuations using this correlation 
function, we combine Q with ( 19 ) and ( 27 1. We then use 
0, (|TTJ>, and (|8]| to calculate v 2 {2} and u 2 {4} in FigjTJ 
where the flow coefficient relative to the reaction plane 
is (v n ) = J pi(p) cosn(</> — $> RP )dp. It is significant that 
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a\ oc 7Z; this is also true for the long range contribution 
to other fluctuation quantities pQ. 

The influence of event shapes and eccentricities en- 
ters (27) in two ways. First, the correlation function 
c(xi, x 2 ) includes a probability distribution of sources 



( 20 1 that implicitly accounts for event-wise variation of 



the volume and eccentricity. Integration over spatial co- 
ordinates effectively represents an average over all possi- 
ble event shapes. Second, the azimuthal distribution of 
particles from each source is modified by flow. A pair 
emerging from a source that experiences a greater push 
will have a narrower opening angle. The magnitude of 
the push follows (25) and depends on position. 



It is important to note that any radial expansion will 
induce correlations contributing to all orders of v n . Par- 
ticular choices of the flow velocity affect their relative 
contributions, i.e. an elliptical flow velocity profile in- 
duces a large i>2 and a triangular flow velocity profile 
would induce a large V3. Equation ( 25 ) is chosen to mimic 
the average (elliptical) flow profile. 

The long range behavior introduced by Glasma corre- 
lations provides key centrality and collision energy de- 
pendence on the saturation scale Q 2 S through 1Z in (23 1. 
Although, the blast wave parameters T, v s , and e also 
provide significant centrality dependence, their change 
with collision energy is minimal. The average v s and 
freeze out temperature in 62 GeV and 200 GeV Au+Au 
collisions are the same as those used in [TJ [FJ |U] and are 
based on an analysis in [33]. At 2.76 TeV the velocity 
is scaled up from the 200 GeV values by 6% and the 
temperature is scaled up by 7% as presented in pQ. 

In Fig(2]we show comparisons of our calculated i>2{2} 
(top panel) and t>2{4} (bottom panel) to measured data 
[H 31 [S] and find the agreement concerning change in col- 
lision energy is quite good. Similarly in Figj3]we calcu- 
late V4{2} and ^4{4}, represented as the solid and dashed 
lines, respectively. Comparisons to u n {4} in both figures 
show how the eccentricity e affects the results. Recall 
that e is tuned to fit in Au+Au 200 GeV collisions. 
The effect of eccentricity is to reduce the even harmonics 
as the collision area becomes circular and e — > 0. The 
non-zero v n {2} in central collisions are due to fluctua- 



tions ( 14 ) , which in our case are due to Glasma correla- 



tions (27). The effect of the change in blast wave param- 



eters with collision energy is apparent from the separa- 
tion of the solid lines in the bottom panel of Figj2] and 
the dashed lines in Fig. [3] The importance of the Glasma 
contribution is similarly evident by the change is separa- 
tion between w„{4} and v n {2} with collision energy. The 
increase in growth of flow fluctuations ([8| is coupled to 
the growth in the saturation scale. 

STAR has also measured vi fluctuations in the form of 



(28) 



resembling the so-called "coefficient of variation" defined 
as the standard deviation divided by the mean. Care 
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} 

0.1 
0.08 
0.06 
0.04 
0.02 




v 2 {4} 
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T ALICE Pb-Pb 2.76 TeV v 2 {2) 

• STAR Au-Au 200 GeV v 2 {2} 

• STAR Au-Au 62.4 GeV v 2 {2} 
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V ALICE Pb-Pb 2.76 TeV v 2 {4} 
□ STAR Au-Au 200 GeV v 2 {4) 
O STAR Au-Au 62.4 GeV v 2 {4) 
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FIG. 2. Measured v 2 {2} (top) from STAR 2 and ALICE 
[H [5] compared to calculations usin g (pi and (|8[) - Same for 
V2{4} (bottom) computed from |T|, ( |11[ ), and ( |24[ ) 



,0.02 




300 350 400 
N part 



FIG. 3. Same calculation as in Fig[2]but for the fourth order 
flow coefficients ?m{2} (solid) and U4{4} (dashed). Data is 
from the ALICE coflaboration [1]. 



should be taken here since the definition of c^, Eq.(14|, 
is not strictly the variance. In Fig|4]we compare our cal- 
culation of (28) to measurement. Calculations at RHIC 



energies seem to agree reasonably well and we include the 
calculation for Pb+Pb 2.76 GeV as a prediction. 

Observe that local correlations that contribute to a n 
lead to fluctuations of all harmonic orders in n including 
v 3 . As a result, Glasma correlations generate a measur- 
able i'3{2} even in the absence of ^{4} or any triangular 
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FIG. 4. Flow fluctuations in terms of the coefficient of varia- 



tion for elliptic flow (28 1 compared to STAR data [2]. 



FIG. 6. The triangular flow coefficient, 113 {2} compared to 
STAR and ALICE data [3j[4]. 



VI. THE RIDGE 
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FIG. 5. Triangular flow fluctuations, a V3 , calculated from 
Q and compared to data computed from t>3{2} and i>3{4} 
measurements by the STAR and ALICE collaborations (3J E] - 



flow. In Fig(5]we show ^3 fluctuations from Glasma. We 
further emphasize that the energy dependence in Fig[5]is 
in good accord with data, supporting the Glasma scaling 
with Q 2 S . The shape with centrality reflects contributions 
not only from Glasma, but also our parameterizations of 
the average reaction plane eccentricity e and our choice 
the flux tube distribution p FT . 



To use ([8| and (14) to calculate vs{2}, we must come 
to grips with the fact that our blast wave parametrization 
assumes i>3{4} = 0. While this seems to be the case for 
the STAR measurements, ALICE has measured a non- 
zero «3{4} for Pb+Pb collisions at 2.76 TeV. To correct 
for this possible discrepancy, we can provide an ad-hoc 
parameterization of «3{4} and use ((8| and (14) to calcu- 
late «3{2}. Agreement shown in Figfej is reasonable. As 
a preferable alternative, we compare our calculated a n to 
fluctuations extracted from ALICE and STAR measure- 
ments in Fig [5] 



Early analyses of the soft or untriggered ridge, includ- 
ing our own, focus on the idea that the ridge is com- 
posed of correlations in excess of momentum conserva- 
tion and elliptic flow. The novel observation of Alver and 
Roland that correlations resulting from geometrical fluc- 
tuations result in odd flow harmonics such as triangular 
flow, V3, suggests that v$, together with flow harmonics of 
all higher orders, explains this excess [55]. In this paper 
we have studied how angular correlations tied to com- 
mon production points influence two-particle correlation 
measurements of flow harmonics v n {2} through flow fluc- 
tuations cr^. Most notably, as discussed in SecJvJ we find 
significant contributions to «3{2} even in the absence of 
triangular flow. In fact, the same correlations contribut- 
ing to flow fluctuations were initially proposed as expla- 
nations of the ridge [7HH UMT7\ HO [50], although their 
role in reproducing the full A<j) azimuthal structure was 
not understood. 

Flow and its fluctuations both contribute to the two- 
particle correlation landscape. The pair distribution is 

p 2 (A0) =Prcf ^l + 2^( t ;„) 2 cosnA^ +r(A0), (29) 

where /? re f is the experimental mixed-event background. 
Observe that the Fourier coefficients of ( 29 1 reproduce 
([7]). We take p rc f = ^ / Pi J Pi where the overbar indi- 
cates an event plane average, to mimick the experimental 
mixed event technique for constructing p rc f. The quan- 
tity Ap(Acj)) = p2(A.(j)) — p rc f characterizes the angular 
correlations. STAR measures the ratio 



Ap(A^) 




cos nAcj) - 



r(A4>) 

Prof 



(30) 

where the (v n ) terms follow from ([lj. 

STAR performs a multi-parameter fit to this distribu- 
tion to separate the ridge peak from the elliptic flow con- 
tributions (along with momentum conservation and HBT 



3 I t ALICE Pb-Pb 2.76 TeV Preliminary 
^p r C ■ STAR Au-Au 200 GeV 




of (fsl) and (31 1, one can write 



6 v=2*l\l bin /N part 



FIG. 7. Near side ridge amplitude calculation from Glasma 
source correlations. Experimental data is from (STAR) |5H - 
[53] and (ALICE) jSJ 



contributions , which we omit for clarity). They report a 
flow-subtracted ridge amplitude (Ap(A</>, A77)/ \ffhtf)\FS 
on the near side, centered at A<fi = Ar/ = [5TH53] , We 
use ( [23| and (27) to calculate the flow-subtracted ridge 
amplitude 



Ap 



dN 

= n d ^F{Acj>), 

fs d y 



(31) 



where F(A(f>) oc r(A0) is the angular correlation function 
normalized so that J F(A<fi)dA(f> = 1. The energy and 
centrality independent scale constant k in ( 23 ) is fixed 
by Au-Au 200 GeV data as in [HUH]- 

The blast wave parameters have little energy depen- 
dence and the Glasma factor ( 23 ) allows for strong agree- 



ment with the 62 GeV data without adjustment of k. Ad- 



ditionally, the ALICE collaboration has measured (31) 
in Pb+Pb collisions at y/s —2.76 TeV [54], providing 
a further test of the CGC-Glasma energy dependence. 
We confront the ALICE data in the same way, without 
adjusting n, however the ALICE [54] measurement pro- 
cedure differs slightly than that from STAR [53 . This 
difference in fitting procedures is mostly insignificant, es- 
pecially in central collisions, but not necessarily so in pe- 
ripheral collisions to which we normalize our calculation 
(at 200 GeV). To adjust for this possible effect, we mul- 
tiply (31 ) by an additional scale factor of ~ 1.5 to return 
agreement in peripheral collisions; see note [55]. In Fig{7] 
we compare calculated ridge amplitudes from (31) with 



measured data for 200 and 62.4 GeV Au+Au collisions 
from STAR as well as 2.76 TeV Pb+Pb collisions from 
ALICE. In this paper we only compare to measurements 
from which elliptic flow has been subtracted, because we 
find that fluctuations dominate w 3 at STAR. 

The characterization of the ridge in terms of Fourier 
coefficients can prove a valuable tool for analysis of such 
discrepancies as well as relate the magnitude its effect 
with respect to other phenomena. Viewing ( 14 ) in terms 



dN 
2 ^ 



A P {Acp) 



cos(nA0) dAcj), 



(32) 



A similar equation was pointed out by Sorensen et. al. 
[551 IS"?] , and discussed in terms of event plane eccentricity 
fluctuations. 



We stress, that (32) is a general result that is model 
independent. Direct measurement of this quantity can 
isolate the effects of flow fluctuations on the ridge by com- 
paring the results to independent flow fluctuation mea- 
surements. Furthermore, flow fluctuation measurements 
can be an effective way to quantify the ridge. Devia- 
tions from (32) at large Arj. can indicate the degree at 



which other phenomena contribute to long range corre- 
lation measurements. At smaller Arj shorter range phe- 
nomenon such as jets or diffusion also come into play. 
Combined flow fluctuation and ridge studies analyzed 
with increasing lower p t limits as done in |58j can reveal 
the emergence of jets and their influence. For example, 
one could test contributions to v n {2} from correlations 
induced by jet quenching as suggested in [31151?]. 



VII. FACTORIZATION 

We now ask whether the azimuthal dependence of the 
pair distribution is entirely determined by the flow co- 
efficients. This is somewhat of a circular question, be- 
cause flow coefficients are themselves derived from cor- 
relation analyses. A meaningful answer depends on how 
the flow coefficients are defined. Experimentalists expand 
the transverse-momentum-dependent pair distribution as 



p 2 oc 1 +2y^V nA (p tl ,p t2 ) cos(nA^). 



(33) 



The term "factorization" is often applied when one can 
express the Fourier coefficients V n A{pti,Pt2) as a product 
of flow coefficients v n (pti)- Authors often cite factoriza- 
tion as a signature of collective flow [20 -22, 60]; see also 
[2"5] for a theoretical perspective. Problems stem from 
the fact that V u a = v n {2} 2 as an exact consequence of 
the definition Q. In other words, V n A always factorizes 
as w„{2} x v n {2}, regardless of the source of anisotropy. 

To define factorization in a useful way, one must com- 
pare V n A to the coefficients (v n (p t i)) defined by the reac- 
tion plane as in (JlJ. Proxies like u„{4} can also be used. 
To see why this is true, we define r„ to be the Fourier 
coefficient of the two particle correlation function r(A</>) 
in We then compare Q to ([29]) to find 



V nA (pti,Pt2) = {v n {pti)){v n {pt 2 )) + r "^ tl,Pt2 ) (34) 

Pmi(Ptl,Pt2) 

where, for simplicity, we keep only the leading order in 
the number of particles in the relevant p t ranges. The 
first term in ( 34 ) measures the effect of global geometric 
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anisotropy, while the second term is due to local fluc- 
tuations. Integrating over momenta and using we 
indeed find 



Vn 



2a 2 n = v n {2} 2 



(35) 



as noted above. However, (34 1 shows that the Fourier 
coefficients of the correlation function r n (p t i,pt2) break 
factorization in terms of (v n (pt2)}- 

We see that hydrodynamic flow does not generally im- 
ply factorization, as follows from (27) and (34). Specifi- 



cally, ( 32 ) shows that flow fluctuations violate factoriza- 
tion. In addition, jets, HBT, and resonances contribute 
to r n along with the long range correlations we have con- 
sidered here. 

Flow coefficients at the LHC are found to factorize for 
pairs of high p t particles, but not when high and low p t 
particles are mixed (60] - This may seem puzzling if you 
view factorization as a hydrodynamic signature. To un- 
derstand this result, observe that the contribution r n /p Te { 



in ( 34 ) can be large compared to anisotropy only when 



a substantial fraction of pairs in the selected momentum 
range are correlated. RHIC energy measurements show 
that the overall magnitude of the correlation function 
drops precipitously relative to mixed-event pairs when a 
lower p t cutoff is increased above ~ 1 GeV [58] . Our cal- 
culations in [9 , which include both flow and jet- medium 
interactions, agree with that trend; see, e.g., fig. 9 of 
[jjj- We therefore do not expect significant violations of 
factorization from the components r n due to long range 
correlations from flow or jets above p t i > 1 — 2 GeV. An 
analogous trend at higher beam energy may explain the 
LHC results. 



VIII. SUMMARY 

In this paper we study the connection between flow 
fluctuations and initial state correlations. Flow fluctua- 



tions are defined by ( 14 ) using two and four-particle cor- 



relation measurements of the harmonic flow coefficients. 
Section [H] provides a general discussion of the calculation 
of flow observables from multiparticle momentum space 
distributions. Following the cumulant expansion method 
discussed in [TBI US], multiparticle distributions can be 
expanded into factorized and correlated parts, e.g. ^ 
and (10); correlated parts contribute directly to the fluc- 



tuation observable a n . 

We study the contribution of local long range correla- 
tions to <r n based on parton production at common trans- 
verse positions. Section [V] presents calculations of a n in 
our Glasma flux tube model. We mention that our model 
incorporates much of the same physics as used by event- 
by-event hydro in [XT] . Those authors obtain the 
initial state in individual events by sampling a probability 
distribution analogous to our p FT . In all cases, Cooper 
Fry freeze out is used. We use a blast wave constrained 
by data to approximate the velocity distribution and the 
freeze-out surface. The key differences are that 1) our 



correlated regions are point-like in the transverse plane, 
while theirs may be larger and 2) they evolve their ini- 
tial distributions using deterministic hydrodynamics in 
each event. We expect that they will eventually achieve 
a better description of the v n at each energy. Our com- 
plementary aim is to study the big picture by exploring 
the energy range with a variety of variables. 

The Glasma formulation provides key collision sys- 
tem, energy, and centrality dependences resulting in rea- 
sonable agreement with experimental measurements of 
v n {2} from 62.4 GeV Au+Au to 2.76 TeV Pb+Pb colli- 
sions for n=2, 3, and 4. The presence of a n provides two 
key results: non-zero values of even harmonics in central 
(circular) collisions, and the existence of ^{2} without 
triangular flow. 

In Sec |VI| we turn to discuss two particle correlations 
and, in particular, the ridge. We see that flow fluctua- 
tions and the ridge are facets of the same phenomenon, 
related by a Fourier transform ( 32 ) . The computed peak 
ridge amplitude in Fig. [7] is in reasonable agreement with 
calculations from 62.4 GeV Au+Au to 2.76 TeV Pb+Pb, 
consistent with our flow results. Further investigation is 
needed to understand the full correlation landscape. In 
particular, we did not discuss the n = 1 harmonic, in 
which momentum conservation plays an important role. 
Characterization of the ridge in terms of flow fluctuations 
(i.e., harmonics) combined with earlier fit procedures as 
in Ref. [S3] will surely prove a useful tool. 

We discuss the question of whether the Fourier coef- 
ficients of the correlation function factorize into a prod- 
uct of flow coefficients in Sec |VII| In general, irreducible 
two-particle correlations from any mechanism violate this 
factorization [2"UH2"2"] . In particular, the long range corre- 
lations that we compute violate factorization at low pt, 
as would jets, resonance decays, and HBT effects. We ex- 
pect any such correlations to become negligible for pairs 
above p t ~ 1 — 2 GeV, due to the rapid decrease of the 
number of correlated pairs relative to the mixed-event 
background that is observed experimentally [58] . 

Finally, we observe that early-time fluctuations have 
broad implications beyond the flow fluctuations and az- 
imuthal correlations studied here. In Ref. [TJ we stud- 
ied the impact of these correlations on multiplicity and 
Pt fluctuations. We argued in Sec. [Hi] that p t and flow 
fluctuations are intimately related because they are both 
driven by local hydrodynamic fluctuations. While the 
average flow coefficients are primarily determined by the 
global event shape, the fluctuations of these coefficients 
have more in common with other fluctuation observables. 
Importantly, the overall magnitude of the contribution 
from long range correlations to all of these fluctuations - 
flow, multiplicity, and pt - is set by a single scale factor 1Z. 
In Glasma theory, the dependence of 1Z on energy, cen- 
trality, and projectile mass are fixed its variation with the 
saturation scale Q s . Results here and in Ref. [JJ provide 
a survey of the ridge, and multiplicity, momentum and 
flow fluctuations that reveals a common energy and cen- 
trality dependence that we attribute to the production 
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mechanism. Glasma calculations are consistent with this 
dependence. 
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Appendix A: Four-Particle Flow Coefficients 

In this appendix we calculate corrections to the four- 
particle flow coefficient measurement v n {4} from two- 
particle correlation sources. Trying to keep the notation 
as general as possible we follow [18 using the cumulant 
expansion method to write 

{e M*-M) = I P2(Pi,P2)e"^-^d Pl dp 2 (Ai) 



/ p2{Pi,P2)dpidp 2 
(N) 2 (v n ) 2 J r( Pl ,p 2 )e mA *d Pl dp 2 



(N(N-l)) 



(N(N-l)) 



(A2) 



where v n {2} 2 = { e m ^~^). Notice that since sin(nA0) 



is odd, the second term in ( A2 1 is equivalent to ^ when 
the pair distribution is symmetric about the normal to 



the reaction plane. Based on the definit ion (21 1 the fac- 
tor (N) 2 /(N(N~ 1)) in the first term of jA2| ) contributes 
a correction of 1/(1 + 72.) w 1 as long as the multiplicity, 
N, is large. For example, in central 200 GeV Au+Au 
collision we calculate 1Z ss 0.003. This number will de- 
crease with increasing collision energy but will increase 
as collisions become more peripheral. This behavior is 
discussed in pQ. 



Our goal is to calculate the four-particle flow coefficient 
«„{4} 4 = 2(e mWl - <fe) ) 2 - (e^^+^-^-W), (A3) 



where the final term is calculated analogously to ( Al ) but 



from the four-particle cumulant expansion (10). Keeping 



only contributions from two-particle correlations we have 

(N(N - 1)(N - 2){N - 3))( e ™(^+^-03-04)^ = ( A4a ) 
= (N) 4 (v n ) 4 (A4b) 

+ (N) 2 (v n ) 2 ■ 2Re y r( Pl ,p 2 )e l2 "^-^p)d Pl dp 2 | 

(A4c) 

+ 4(A) 2 (i; n ) 2 J r(p 1 ,p 2 )e mA Up 1 dp 2 (A4d) 
fr(pi,p 2 )e l2n (*-^5 d p ld p 2 (A4e) 

J r( Pl ,p 2 )e" lA ^ Pl dp 2 . (A4f) 

where = <j>i — <\> 2 and $ = (4>i + <f> 2 )/2 are the relative 
and average coordinates. To reduce the equations further 
it is necessary to make the approximation that (N (N — 
1)(N -2){N- 3)) « (N(N - l)) 2 then the terms ( |A4d[ ) 
and ( A4f[ ) will cancel with corresponding terms emerging 
from the twice the square of ( A2 ) , leaving the corrections 



VnW = K> 4 - K> 2 • 2i?e {S 2 } - |S 2 | 2 , (A5) 
where 



Jr(p 1 ,p 2 ) e ' 2 "(^-^)dp 1 dp 2 
(N(N - 1)) 



(A6) 



which is equivalent to (13 1. Notice that these corrections 



are effectively of order 2n and depend on the reaction 
plane. As discussed in the text, the corrections in (A5) 
maximally modify (11) by 1.2% in our model. 
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